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ABSTRACT 


The  motion  of  a  projectile  In  a  gun  barrel  la  prescribed  so  that 
the  motion  of  the  gas  In  the  part  of  the  barrel  adjacent  to  the  pro¬ 
jectile  Is  given  by  a  simple  analytic  function.  Analytic  methods  In¬ 
volving  the  hodograph  transformation  are  employed  to  determine  the 
motion  of  the  gas  In  the  remaining  portion  of  the  barrel.  The  pro¬ 
cedures  are  carried  out  explicitly  in  the  case  y  ■  — 


FOREWORD 


The  investigations  of  this  report  were  performed  in  the  Applied 
Mathematics  Section  of  the  Mathematics  Research  Group  as  a  project 
supported  under  WEPTASK  No.  WR-4-0046.  The  investigations  were  under¬ 
taken  in  connection  with  a  problem  proposed  by  Dr.  Arnold  E.  Seigel  of 
the  U.  S.  Naval  Ordnance  Laboratory.  The  date  of  completion  was 


I.  INTRODUCTION 


In  the  problem  to  be  considered, 
a  projectile  P  moves  through  a  gun 
barrel  with  a  prescribed  motion, 
describing  the  path  OPQ  In  the 
position  vs.  time  plane.  The  motion 
of  the  gas  is  prescribed  behind  the 
projectile  and  at  the  base  of  the 
barrel  such  that  the  gas  equations 
have  a  simple  analytic  solution,  as  long  as  the  gas  velocity  remains  below 
the  local  sound  speed  at  the  base  of  the  barrel  (segment  OM  in  the  diagram). 
At  point  M,  the  gas  velocity  equals  the  local  sound  speed;  it  is  prescribed 
that  thereafter  they  remain  equal  (but  not  constant)  at  the  base  of  the 
barrel  (segment  MR).  The  Initial  analytic  solution  is  regarded  as  holding 
throughout  the  region  QOMN  bounded  by  the  projectile  path  and  by  the  C  - 
characteristic  MN  through  M. 

The  problem  then  is  to  determine  the  motion  of  the  gas  in  the  remaining 
region  NMR.  An  attempt  is  made  in  the  following  to  obtain  an  analytic  so¬ 
lution  in  this  region;  the  details  are  carried  out  in  the  particular  case 
y  =  5/3.  It  remains  to  be  seen  whether  this  solution  is  of  unrestricted 
or  of  limited  validity  in  this  region,  due  to  the  occurrence  of  shocks 
or  other  discontinuities. 

H .  DESCRIPTION  OF  PROBLEM 

The  motion  of  the  projectile,  in  normalized  form,  is 


1 


for  position  and  velocity  respectively.  (The  original  unnormalized  problem, 
together  with  conversion  factors,  is  stated  in  the  appendix.)  A  polytropic 


gas  lies  behind  the  projectile;  along  the  projectile  path  its  pressure  and 


sound  speed  remain  at  their  Initial  values: 
(3)  ^  =  /*.  -  /  , 


In  normalized  form. 

At  the  base  of  the  barrel,  x  ■  0,  the  gas  velocity  Is  again 
(5)  t*  --  ^ 

as  long  as  this  remains  less  than  the  local  sound  speed. 

At  time  t  ■  (point  M) ,  the  velocity  equals  the  sound  speed;  specifically 


we  then  have 


(7)  ^  = r 

Is  specified  at  x  •  0  (along  MR).  (Note  that  this  problem  Is  meaningful  for 
a  polytropic  gas  for  1  <  y  <  3  only;  this  assumption  Is  made  throughout.) 

The  continuity  and  Euler  gas  equations  are  respectively 


(8)  <7^  +  * 

and 

(9)  y-  U 
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Here,  for  a  polytropic  gas, 


(10)  * 

is  the  Riemann  quantity,  such  that 

(11)  M  y-  cn  T 

and 

(12)  M  7 

are  invariant  along  and  C~  -character istics  respectively. 

With  our  assumptions,  we  may  take 

(13)  U  -  ^ 

throughout  the  region  QOMN.  The  gas  equations  then  reduce  to 

(14)  crif  v-  ^  — *  * 

and 

U5)  /  v-  2-^ 

these  have  the  general  solutions 

(16)  ^  f  fX - ZT~ / 

and 

(17)  S'*  + 

(f(z)  and  g(z)  arbitrary  functions)  respectively.  The  resulting  solution 
in  QOMN  is 

da)  <r-  f  */ =  /? -  ■ 

«-2  JZ  j 

(At  the  projectile,  x  =  ,  we  have  «-  »  /  .  Hence  by  (17), 

•  111  en  by  (17>  again, 

,  //  ,  ^  TjEv  E  jT  ,  which  agrees  with  (16).) 

<^=  y  f--/  (  *  /  r-' 
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The  C“ -characteristic  MN  through  M  must  satisfy 


(19)  //- 


”  ***  —r  /V  “  ? — / 


this  yields 


<‘°>  '  *  4^7  -  *?/>-  s=-v®7 

as  the  equatior  for  MN. 

The  problem  now  becomes  that  of  determining  u(x,t)  and  ^*(x,t)  in  the 
region  NMR,  such  that  they  satisfy  the  gas  equations  (8)  and  (9)  within 
the  region,  and  such  that  the  initial  conditions 
(21)  //  sr 

hold  along  x  “  0  and  (the  consistent  conditions) 


(22)  U  —  ^ 

hold  along  X  *r 


III.  METHOD  OF  SOLUTION 


The  solution  of  the  problem  described  above  will  be  attempted  by 
means  of  the  hodograph  transformation.  From  (11)  and  (12),  we  obtain 


(23)  //  *■  JT-  7  , 

(24)  O'--  T+  7  , 

in  terms  of  the  characteristic  coordinates  (T,r>»  where  JT  is  a  constant 
along  any  C+ -characteristic  and  is  a  constant  along  any  C“ -character istic . 
In  these  coordinates,  the  gas  equations  become  equivalent  to  the  equa¬ 


tions 

(25)  xr  -  ^ 


and 
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(2  6)  Xy 

Let 


-2 


* 


(27)  A  - 


■*  /  *=7/ S 


upon  differentiation  of  the  above  equations  with  respect  to  and  JT* 
respectively,  and  elimination  of  ,  a  linear  equation  for  ^ 

alone  results: 

(28)  ^ZV 


f? 


Z2/  ^  *?/ 


In  the  ( jT,^r)-plane,  the  C~  and 
C+ -characteristics  are  horizontal 
and  vertical  lines  respectively, 
while  the  region  NMR  now  is 


0  . 


bounded  by  the  straight  lines  I 

MN:  7  ~  and  MR: 


The 


(2  9)  X  =  * 


along  line  MR,  and  (the  consistent  conditions) 

do  r-  ■*  =■  ^£r-  ^ /£££]* 

along  the  C“-characteristic  MN. 

Equation  (28)  for  -**■(/*,  7)  is  one  for  which  the  Riemann  function  is 
known  (Reference  2  !  pp.449-4Q),  thus  enabling  the  solutions  of  normal 
initial-value  problems  to  be  written  in  explicit  form.  The  Riemann  function 
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may  be  written  in  cither  of  the  forma 


(31)  /fj  {*;>-/  -t  7S  -  f  jtrfj  A;  /;  - 

_  ( *  Fr.  \  .  /  .  r*~  Wr-  zS  7 

L  '  '  '  f*  +  7//y*&  S 


where  F(a,b;c;z)  is  the  Gaussian  hypergeoine  tr  ic  function.  (Note  that 

02)  /v^  *s  =  f /-  Jy 

from  which  the  equivalence  of  the  two  forms  follows.)  When  A  Is  an 
integer,  and  thus  in  particular  for  the  important  cases  *■  ^  {A  =*--?/ 
and  ^  *=•  — —  /.  ¥~  f  4  *=■  »  the  Riemann  function  is  rational 

in  its  four  variables.  Specifically,  for  »  Hp  we  find  that 


(33)  y ;  f 3**- 7/ * 


When  >4  is  an  integer,  we  may  also  write  the  general  solution  of  (28) 
explicitely  in  the  form  (Reference  [ 1 1 ,  p.  90) 


(34) 


JC  f 


■zA-'  f-ffr/  7 
tfr'i’J-'  J 


»7 


rr  /?s  7 

-'Lfr+ryy. 


(There  the  constant  and  the  functions  f(z)  and  g(z)  are  arbitrary.  This 
furnishes  an  alternate  method  of  solution  for  these  cases. 

Our  initial -val ue  problem  is  abnormal  in  that  it  is  partly  characteristic 
and  partly  non-characteristic,  and  in  that  only  X  (JT>^)  is  known  along  the 
non-characteristic  part  of  the  boundary.  iliis  complicates  the  method  ot 
solution  if  the  general  Riemann  method  is  to  be  used  to  solve  (28)  directly. 
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IV.  THE  SOLUTION  FOR  THE  CASE 


1 


5/3 


The  solution  has  been  carried  out  explicitely  tor  the  case 
/”  ■  yr  fA  *■  »  the  simplest  non-trivial  case.  The  methods  used 

here  could  be  extended  immediately  to  the  other  cases  with  integral  A  . 

The  Riemann  function  method  in  principle  is  extendable  to  the  case  with  A 
arbitrary,  but  some  other  method  for  solving  the  integro-dif ferential 
equation  (75)  for  T(^),  i.e.  t(/*,  fr  )  along  MR,  would  appear  to  bf  re- 
quired  then. 

The  solution  obtained  here  is  in  the  hodograph  plane  and  is  valid 
throughout  all  of  NMR  there.  However  the  range  of  its  validity  in  the 
original  plane  remains  to  be  determined;  this  depends  on  the  mapping  into 
the  original  plane  remaining  one-to-one,  i.e.  on  the  characteristics  of  the 
same  family  not  intersecting  within  NMR  in  the  (x,t)-plane.  This  can  readily 
be  investigated  numerically  by  using  the  solutions  derived  here;  a  short 
computer  program  for  this  purpose  is  now  being  written.  (The  line  MR  in 
the  original  plane  is  known  to  be  singular,  however.) 

The  solution  was  first  obtained  by  using  the  general  solution  (34). 

As  this  method  is  limited  to  Integral  A  ,  an  attempt  was  then  made  to  use 
the  Riemann  function  method.  With  the  latter  method,  the  problem  separated 
into  two  steps;  the  first  step  led  to  an  integro-dif ferential  equation 
essentially  for  the  missing  initial  value  t(jT, jr)  along  MR.  The  second 
step  is  then  the  ordinary  non-characteristic  Riemann  solution,  with  the 
above  Initial  value;  this  is  straightforward  and  is  not  carried  out  here 
in  detail  (see  Reference  [2],  pp.  449-453). 
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A.  General  Solution  Mathod 

For  ^  ,  equation*  (25),  (26)  and  (28)  bacon* 


os)  xr  -  ? 

(36)  Xf  -  f  fax-  ?S 
and 

(37)  y-  t +■**/-  *■ 

The  initial  condition*  become 

(38)  X  ~  & 

along  line  MR:  jf  ^  ^ 

and 


(39) 


along  NM:  ^ 


S7 

x  y 


-r 


* 


The  point  M  ha*  coordinate* 


(40) 


j  ^ 


here  the  initial  values  are 


(41)  AT  *  0  y  ^  f 


From  (34),  the  general  aolution  for  t(JT,  jr)  i* 

(«)^4T 

the  constant  of  (34)  turn*  out  to  be  superfluous. 

Using  (35)  and  (36),  we  find  for  x(JT,  p)  the  expression* 

(43)  X  —  *SXT*S*  *■  ? 

=  c'+  f-ffr-x?;  XYZ  f$%*l 
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and 

(44)  X  —  Sx7*/7  *■ 


-  C 


y- 


7 

/>V  V-'V 


(Here  C7  and  C"  are  arbitrary  functions  of  and  respectively.) 

The  two  expressions  must  be  the  same,  hence  C'  =  C"  =  C,  a  constant. 
Working  out  the  details,  we  obtain  finally 

(45)  jc/Jv  v  =  fc(rT+?;J+ 


y-  />jr-  Y?/J-  ^ 


(The  constant  C  will  later  turn  out  to  be  superfluous.) 


On  MR:  ]T  we  have  from  (38)  the  result 


and  so  the  differential  equation  for  g  (?): 

(47)  ^  v?/  -  jjr  *‘''7/  «  -  +  jy  T^^r/. 

The  general  solution  of  this  equation  is 

(48)  sW-  *7  +  “Sr-  “ft**/*. 

/r 

On  NM:  ^  *  ,  we  have 

(49)  ~  r~  ~f~  * 
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and 


(so  xfr,  *fS-  rfr-  & 


-£Js  YfJsfr- 


This  pair  of  linear  algebraic  equations  for  f(.jT)  and  f,(jT)  may  be 
solved  to  give 


(51)  frj, f/*  Ygy- 
and 

(52)  t cfr/-  fir* ■£/£?*: 


this  solution  may  be  verified  by  differentiating  (51), 


The  constants 
may  set 
(53)  n  * 


and  Y-f/  prove  to  be  superfluous,  and  we 


Then,  substituting  (51)  with  r  into  (48)  gives 

(54)  /r 7  '/j -  -£  a  sy yr 1 


J*  **  —  *7* _ & 

7  ~  jf-  v 


and  differentiating  gives 


(55) 


Y?Y  -  y-  J/>-  IT  cfs?- ST/+  7r  ?  S-  # 


ST 


(56)  A k 


/T 

Setting  ^  and  using  (53)  in  either  of  the  above  yields 

7  '/ 


JTJ*0  r  S 
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C  proves  to  be  superfluous  and  may  be  set  equal  to  zero. 


Thus  we  obtain: 


(57)  f/ry-  Jrfr- 


^  __ '  r-'-i’  _  */^r~  f 

j  J  -  f*  •*  y/  V 


(58)  -  s7 


(59) 


'/?/= 


77  ?  ~  3^?  '  77 


(60) 


_/  */ZF~  77  ~  ^  /y  M  ^  M  _ 


31-Z0 


From  (42)  and  (45)  we  now  obtain  the  solutions  within  NMR: 


(61) 


fr+-  vS'i 


**/*■*+ r 


v  *  •- 

//  / 


v-  -?y>  frJ 


#■?*/-  77  ^Jy 


(62)  -y/ fy  7/~  jjT^frf-Trr  7  ^ 

B  .  Rlemann  Function  Method 


yr_ 


-Vj 


J-r'+  ry-j?7J'fr-dzJ. 


Rlemann's  solution  for 


the  noncharacteristic  Initial- 


value  problem  gives  the  solu¬ 
tion  at  a  point  N  in 

terms  of  the  values  of  t  and 
its  derivatives  along  MR.  Our 
problem  is  an  inverse  one,  in 
that  the  value  of  t(jTtp')  at 


V*1? 


U-2V 


N  is  known  while  that  along  MR  Is  partly  unknown.  The  result  is  that  the 


solution  formula  will  constitute  an  integro-dif f erential  equation  for  the 


function  t  along  MR 
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Once  t  has  thus  been  determined,  the  Riemann  solution  may  again  be 
employed,  this  time  directly,  to  obtain  the  solution  for  arbitrary  points 
inside  NMR. 

The  symmetric  form  of  the  Riemann  solution  for  our  problem  is 

(63)  rr/'G  (■  -j*1'  jy 

y  A  fsrs  7/  frts  v/J  y 

where  I(jO  is  the  line  integral  along  MR: 

(64)  -zy>y  /iyyZ  ~ 

where  Rg  is  the  Riemann  function 

(65)  Ajt  *  ■*  7s 

and  ^fi/4  y  ^  for  short. 

(See  Reference  [”2],  page  453,  equation  (4),  where  however  the  term  ^  fiy 
should  be  just  a  A  ,) 

The  line  integral  may  be  evaluated  by  setting 

(66)  «  =■  -  7-  y  *  7~  *  ?  . 

This  gives 

(67)  -T/*/  fcr'  rs -*z  f-*r' r/-  Uis&x-rsJ 

+[/p  r,-*?,  r?*t, 
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The  function  t(2  7",  T')  along  MR  is  essentially  a  function  of  one 
variable;  hence  set 

(68)  7~ f  r/. 

Along  MR  we  ;  re  given  that  (38) 

(69)  X  / 

hence  also 

(70)  X  fa 7^  rS  =  X  X*  far,  rX t-J  m  . 

From  the  gas  equations,  (35)  givea 

(71)  Xu  far,  X/  m  0, 

and  (36)  and  (70)  together  then  give 

(72)  V  ^  ^ 

Thus,  since  7"  is  not  identically  zero, 

(73)  ~  ' 


Therefore , 
(74)  TV*-; 


XfJ7>  f^ry7-y-r  r-/  * 


Thus  the  knowledge  of  x(2T,r)  enables  us  to  replace  the  derivatives  of 
t  along  MR  by  the  single  derivative  T/(7");  (63)  then  becomes  finally 


(75) 


1 


-  A 


f  /&  T-;Jr,  7-YrS^fyL  A  ry  -?*  #/ 

fir,  r;  FJ  +  i  &r'  r'  ■**'  7~' YrsJSr- . 


This  is  the  desired  integro-dif f erential  equation  for  T(jr). 
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The  Riemann  function  and  its  derivatives  for  A  ■  2  are  given  ex- 


pl icitely  by  (33)  : 
(76)  ^ 


and 

(78)  /tsrfrzr;**,  & 


'  fi£zp£r7’*+*/!*r‘ 

r*T~ 


Also,  we  have  from  (39)  the  values 
(79)  7y  *■  ^7  "  ^ 


and 

(80)  7-f-f/  -  fTs 


Equation  (75)  can  then  be  written  in  the  form 


*rj  t/r rj  7~ fr-/?  T- 


The  term  in  T/(7")  can  be  eliminated  under  the  integral  sign  by  an 
integration  by  parts;  this  yields 

ftp 

(82)  f.  / <"r  7-V»^>r/r//VVr«/'f?i.  ^ 7/^y 

* 

-  -  ^£l/rr*  +  ir/r+sr 7]7Yt-/ Sr- . 
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Substituting  Into  (81)  gives  the  integral  equation 


As  the  coefficient  of  T(T')  under  the  integral  sign  is  a  polynomial 
in  (83)  can  be  reduced  to  a  differential  equation  for  T(^f)  by 
successive  differentiations  with  respect  to  Here  two  differentiations 

are  required,  which  yield  the  Euler  equation 

(84)  7-  W  +  7T^  'W  +  jf'  W  -  TF  . 

As  initial  conditions,  we  have  (80),  and  from  (83)  once  differentiated, 
also 


The  solution  of  (84)  satisfying  (80)  and  (85)  is 

(86)  7~/7/  “  7T?  '  7?  7~'~  7T7  7  f 

This  is  the  required  value  of  t  along  MR;  it  can  be  readily  verified  that 
it  is  identical  with  (61)  with  2  y  substituted  for  JT , 

The  Riemann  solution,  with  (86),  03)  and  (74),  may  now  be  used  in  the 
usual  manner  to  obtain  t (Jr*,  ^r)  at  arbitrary  points  within  NMR. 

In  the  general  case,  where  the  Riemann  function  does  not  reduce  to  a 
rational  function,  the  principal  requirement  for  a  solution  by  this  method 
would  appear  to  be  a  practical  method  for  solving  the  Integral  equation 
corresponding  to  (83). 
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A  short  computer  program  has  been  written  to  evaluate  and  plot  the 
formulas  for  x  and  t  obtained  In  the  preceding  section.  A  plot  where 
the  data  has  been  run  out  to  /  ■  2»^and  ^  ^  follows,  with  steps 

of  A  &  &  nd  A  J  .  No  singularities  aside  from  the  singular 

line  x  ■  0  appear  in  the  plot;  the  C_ -characteristics  are  all  tangent  to 
the  line  at  x  ■  0.  (In  reading  the  plot,  it  should  be  noticed  that  an 
additional  vertical  line  of  points  falls  almost  on  top  of  the  line  x  =  0, 
so  that  the  parabola-shaped  C~ -characteristics  actually  end  one  point 
below  on  x  ■  0  from  their  apparent  end  points.) 

A  plot  with  the  data  run  out  to  JF  ~ 
is  very  similar  in  appearance,  with  no  singularities  aside  from  x  *  0 
apparent.  However,  because  of  the  non-linearity  of  the  problem,  it 
cannot  be  assumed  from  this  alone  that  singularities  inside  the  region 
of  interest  will  not  eventually  appear.  Additional  runs  with  data  run 
out  to  increasingly  greater  distances  are  planned. 


16 


T(XSI.ETA)  vs.  X06I.CTA) 


REFERENCES 


1  R.  Courant  and  K.  0.  Friedricks,  Supersonic  Flow  and  Shock  Waves. 
Interscience  Publishers,  1948,  pp.  79-91,  191-197. 

2  R.  Courant  and  D.  Hilbert,  Methods  of  Mathematical  Physics.  Volume  II, 
"Partial  Differential  Equations".  Intcrscience  Publishers,  1962, 

pp.  429-432,  449-461. 


17 


APPENDIX  A 


APPENDIX  A 


The  Unnormalized  Problem 


The  problem  as  stated  in  the  body  of  the  report  is  in  fully  normalized 
form.  The  original  form  of  the  problem  is  stated  below,  following  which  is 
given  a  set  of  conversion  ratios  between  the  unnormalized  and  the  normalized 
quantities . 

The  motion  of  the  projectile  is  given  as 


(Al)  *  ^  ^7  ^  y 

(A2 )  4/  -  ^  / 


with 

(A3) 


sr 


0 


Here  is  the  initial  gas  pressure,  A  the  cross-section  area  of  the 
barrel,  and  M  the  mass  of  the  projectile,  in  consistent  units. 

Let  aQ  be  the  initial  sound  speed;  at  the  base  of  the  barrel  the 
gas  motion  is  given  by 
(A4)  y 

(A5)  -5=-' 


as  long  as  M  *<  A  ,  the  local  sound  speed.  At  time 

(A6)  ^ 

the  condition 
(A7 )  *  *■  <*■ 

is  attained;  this  condition  holds  at  the  base  of  the  barrel  for  t  >  t^. 


The  C“  -characteristic  NM  bounding  the  prescribed  solution  la  given  by 


the  equation 


(A8)  X 


X*  * 


rZEZ  /  7  ' 
**  *■+  -/  / 


along  thla  characterlatlc  the  condltlona 


(A9)  // 


(A10) 


hold. 


If  f  represents  a  given  unnormalized  quantity,  let  f  denote  its 
normalized  counterpart  (the  quantity  used  elsewhere  in  the  report),  and 
f  the  conversion  factor  connecting  the  two,  so  that 

(aii)  ** 


Define 


(A12)  ^ 


The  conversion  factors  are  now  the  following: 


(A13)  ^ 


(A14)  X 


(A15)  ~  *  —  cr  * 
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